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Abstract 

Zero temperature properties of a dilute weakly interacting d-dimensional Bose gas in a random 
potential are studied. We calculate geometrical and energetic characteristics of the localized state 
of a gas confined in a large box or in a harmonic trap. Different regimes of the localized state are 
found depending on the ratio of two characteristic length scales of the disorder, the Larkin length 
and the disorder correlation length. Repulsing bosons confined in a large box with average density 
n well below a critical value ric are trapped in deep potential wells of extension much smaller than 
distance between them. Tunneling between these wells is exponentially small. The ground state of 
such a gas is a random singlet with no long-range phase correlation For n > nc repulsion between 
particles overcomes the disorder and the gas transits from the localized to a coherent superfluid 
state. The critical density nc is calculated in terms of the disorder parameters and the interaction 
strength. For atoms in traps four different regimes are found, only one of it is superfluid. The 
theory is extended to lower (1 and 2) dimensions. Its quantitative predictions can be checked in 
experiments with ultracold atomic gases and other Bose-systems. 

PACS numbers: 03.75Hh, 03.75.Kk 
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I. INTRODUCTION 



Bose-Einstein condensation (BEC) is one of the key quantum phenomena which shaped 
the physics of the 20th century. In BEC a macroscopic part of the atoms is in the state 
with lowest energy. Well known examples are the superfluid phase of ^He and laser cooled 
atoms in magnetic or optical traps, which are a playground for the demonstration of generic 

aiantum phenomena like the appearance of quantized vortices, the Josephson effect etc. 
.B.a. 

A more complicated situation arises if these phenomena are studied in a random environ- 
ment like in porous media. Experimentally, superfluidity was shown to persist in random 
media like Vycor or aerogel glass {4] which contain pores filling from 30% up to 99% of the 
material. Partial coverage of the pores leads to a reduction of the transition temperature: 
the transition temperature ranks from 2K for full-pore Vycor to about bmK for low cover- 
age films in Vycor. Finally, for very low coverage superfluidity corresponding to less than 
1.5 monolayers the superfluidity vanishes completely even at T = 0. One must presume 
that the helium is strongly localized there. The remaining Helium atoms were considered 
to form an "inert" layer which screens the random potential once the density of the Helium 
atoms is again increased. 

The experimental investigations on BEC in disordered environments concentrated on 
critical behavior of the superfluid density and specific heat at decreasing coverage of the 
porous media More recently, Bose-Einstein condensates have been studied in optical 
traps which superpose a harmonic and a random potential Most of the 

papers concentrate on the possibility of Anderson localization in effectively one-dimensional 
condensates. In contrast to random potentials in crystals, the random speckle potential has 
a typical correlations length of the order of /im which is about of the same order as the 
transverse extension of the atomic cloud. The longitudinal extension is about 10 to 100 
.a.ge. fl. In .os. of .he experiments, the expanding dens.ty prohle afte. switching 
off the harmonic potential was analyzed. For weak random potential Lye et al. 5|] found 
essentially the density profile corresponding to the Thomas-Fermi theory. For larger values 
of the random potential these authors found a stripe-like pattern of about 2yum spacing and 
damped dipole oscillations. For very strong disorder the atoms are localized in the minima 
of the random potential. Chen et al. lOj also studied in situ and found fragmentation 
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of the BEC. In some papers the superposition of a periodic plus a random potential was 
considered to enrich the number of possible phases [y, |7| . Schulte et al. Gj found features 
in the expanding condensate like pronounced fringes and changes in the axial size of the 
central peak which they trace back to the random potential. 

Theoretical investigations of Bose gases or fluids in a disordered environment focused on 
three issues: 

(i) The critical behavior at the normal to superfluid transition in the vicinity of the transi- 



tion temperature Tc(n) 
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I3I . [1^ , where n denotes the boson density. For fixed disorder 



strength Tc{n) decreases with decreasing boson density n and eventually vanishes for n — > ric- 
(ii) The description of zero-temperature quantum phase transition from the superfluid to 
a non-superfluid (glassy) phase taking place at the critical value ric of the boson density 

HQ- 
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(iii) The microscopic description of the superfluid density (and the condensate) in a disor- 
dered environment in the region of comparatively high average density of the gas n ^ Uc 



and T <^ TJn) 
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26|. 



The problem of the Bose gas in the random environment attracted so much attention 
because it connects two central ideas of the condensed matter theory: the Bose condensation 
and localization. Especially the problem of the interplay of the interaction and disorder is 
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281], but only in the case 



intriguing. So far this problem was studied for Fermi gas 
of a weak disorder it could be done rigorously. Therefore the delocalization induced by 
interaction remains beyond the frameworks of the theory, though most of theorists are sure 
that it does exist. In the case of the Bose gas the situation is even more involved since there 
is no Pauli principle which delocalizes the ideal Fermi gas in 3 dimensions. Our study shows 
that the interaction delocalizes the Bose gas at some critical density. 

Since in this paper we restrict ourselves mainly to zero temperature, below we give a brief 
review only of the works mentioned in (ii) and (iii). 

(ii) The quantum phase transition, taking place at zero temperature and n = ric, exhibits 
in general a critical behavior different from that observed at finite Tc{n). This is because in 
the quantum domain static and dynamic critical behavior are intertwined. This quantum 
phase transition was considered by Ma, Halperin and Lee [Q who mapped the problem onto 
a. d + 1 dimensional classical 0^ theory with disorder completely correlated in time. Their 
description gives rise to new critical exponents which were calculated by e-expansion around 
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d = 4: dimensions. This work has been criticized by Fisher et ah IGj. One conclusion of 



the work of Fisher et al. 



161 ] is that the dynamic critical exponent z should be equal to the 



dimension d of the system in all space dimensions. This result has been disputed recently 
171]. Indeed, quantum Monte Carlo studies show z ^ lA for two-dimensional system, in 



discrepancy with the relation z = d 361 ]. 



Some arguments of the work by Fisher et al. 16| rely on the treatment of Giamarchi and 



Schulz 



181 ] who found the superfiuid to insulator transition in one dimension. This transition 



happens at arbitrary disorder by increasing the interaction to a sufficiently large value. The 
fact that the transition into the non-superfluid phase happens at increasing interaction 
seems paradoxical, but it can be explained by the increase of quantum fluctuations of the 
phase. As we will show, for weak interaction there is a second transition resulting from the 
competition between the disorder and the interaction. We believe that namely this transition 
has its counterparts in higher dimension and is observed in experiments with cooled gases, 
(iii) The fate of the inert layer and deeply localized (Lifshitz) states was discussed by Lee 
and Gunn [2^ without speciflc predictions to be checked by experiments. More recently the 
reduction of the superfluid and the condensate density deeply in the superfluid phase was 



considered microscopically by Huang and Meng [21J. These authors used the Bogoliubov 
transformation and found a decrease of the superfluid density and of the condensate if disor- 
der is taken into account. The corrections to both quantities are proportional to — {ric/nf'^'^ 
in our terms. Thus this approach is restricted to weak disorder and n 3> Uc- These results 



have been conflrmed and extended by other authors 22 



23 



24 



25 



261 . All these works 



used the approximation of weak disorder. The extrapolation of their results to the range of 
strong disorder made in 2]| was illegitimate. 



Zhou 



argued 



that the T = transition to the superfluid phase is of the percolation 
considered a BEG in a special potential which is random in one and 



type. Lugan et al. 

parabolic in the two remaining directions. The disorder is bounded from below but not from 
above. Using the numerical solution of the Gross-Pitaevskii equation as well as a mean fleld 
like approach they identifled three different regimes: the Lifshitz glass, the fragmented BEG 
and the non-fragmented BEG. Because of the difference in the models (our potential includes 
randomness in all directions) a direct comparison with our results derived below is not 
possible, except for the one-dimensional case when the harmonic trap potential is negligible. 
Only in this case, and using our characteristics of the random potential, it can be shown that 
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their approach gives results for the fragmented BEC comparable to ours. However we do 
not see a difference between the Lifshitz glass and the fragmented BEC. Moreover, Lugan et 
al. 8| did not obtain important geometrical and energetic characteristics of the fragmented 
state vs. the average density of particles as we did based on a simple physical picture of the 
fragmentation. 

More recent papers considered the dynamic consequences of Anderson localization in 



BEC and do not directly overlap with our work 
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those used in a published work by two of the authors 3J| , Shklovskii 35| independently and 



33|. Usin g ar guments similar to 



practically simultaneously arrived at the same conclusions about the critical value of density 
at which the transition between an insulating and the superfluid phase takes place, but he 
did not give the description of the localized phase. 

The aim of this paper is to develop an alternative approach to the quantum phase transi- 
tion starting from deeply localized state. We deliver a detailed study of the deeply localized 
state of the Bose gas in a random potential. Such a state appears at sufficiently low average 
density. We show, that in this region the Bose-Einstein condensate decays into remote frag- 
ments of small sizes. We give a rather transparent qualitative and quantitative geometrical 
description of the localized state. We show that, at a critical density Uc which we express 
in terms of the disorder characteristics and interaction strength, the increasing tunneling of 
particles between fragments leads to transition from the random singlet state to the coher- 
ent superfluid. The localized state can exist in several different regimes. These regimes as 
well as their quantitative characteristics strongly depend on the ratio of two basic length 
scales characterizing the disorder: the Larkin length C and the correlation length h. We 
extended our consideration to a practically important situation of the Bose gas confined by 
a harmonic trap. In this case the gas forms a cloud which size is determined by the energy 
minimization. We demonstrate that all regimes which appear for the gas in the box appear 
also in the trap driven by the number of particles and other parameters. Our considera- 



tion is based on theory of deeply localized sing! 
potential given in seminal works by 1. Lifshitz 



^e:T)article states in an uncorrelated random 



and Lax 



391, Zittartz and Langer 40[, and Halperin 



411 . Their predecessors Keldysh and Proshko 43|] and Kane 4J] studied electron 



density of state in a smooth random potential ha ving in mind semiconductors. These works 



were extended and detailized in References 



42 



45 



46 



47 



48|, |49|, 150[. In this work we 



propose a modest extension and simplification of their arguments to find some geometrical 
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and physical characteristics of these state necessary for our purposes. On the other hand, 
our theory can be considered as a generahzation of their instanton-type theory incorporating 
the self-consistent field of the interacting particles. 

We consider the ground state of the system at zero temperature. It is not always reach- 
able since the system may be frozen in a metastable state. It is therefore important that 
nevertheless the ground state occurs to be reachable in a sufficiently wide range of param- 
eters. We also demonstrate that properties of the localized states of the Bose gas are rich 
enough and that they are available for the experimental study. 

Our article is organized as follows. In the next section we characterize the random 
potential. In Section HTTl we consider the single-particle deeply localized states in 3 dimensions 
for the uncorrelated random potential obeying the Gaussian distribution. In Section IIVI 
we consider how the dilute Bose gas fills the deep fluctuation potential wells in this case. 
In Section |V] the transition from the random singlet to the coherent superfiuid state is 
considered. Section IVII extends all these results to the case of strongly correlated random 
potential. In Section IVIII we consider the Bose gas subjected to simultaneous action of a 
harmonic trap and random potential. Section IVIIII represents extension of these results to 
lower dimensions 1 and 2. The discussion of our results and conclusions are left for Section 

m 

II. DESCRIPTION OF THE DISORDER 

The disorder will be represented by a random potential U (x) with zero average at each 
point of the space and obeying the Gaussian distribution: 

dW [U (x) , dU (x)] = exp l^-i y U (x) K^^ (x, x') U (x') rfxdx'^ (1) 
X VdetK Yl (Al]/27r)^/^ dU (x) , 

X 

where AQ is the volume of an infinitesimal cell; K (x, x') is the correlation function or 
correlator of the random potential: 

ir(x,x') = (f/(x)f/(x')) (2) 

and (x, x') is the inverse correlator, defined by equation: 

j K (x, x" ) K-^ (x" , x') rfx" = 5 (x - x') . (3) 
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Figure 1: Uncorrelated potential, the characteristic length and energy scales are L and respec- 
tively. There is typically only one bound state in a single potential well. 



The simplest random potential with zero correlation length has a 5— like correlator: 

J^o (x, x') = (t/ (x) V (x')) = (x - x') . (4) 

We will call it uncorrelated random potential. The corresponding distribution function is 
the product of independent distribution functions at each point of space: 



dW^ \U (x) , dU (x)] = exp (^-^ y f/2 (x) rfx^ J] 



27rfi; 



.■2 



dU fx) 



(5) 



More general Gaussian random potential has a finite correlation length h. The simplest real- 
ization of the random Gaussian potential with the finite correlation length is the Ornstein- 
Zernike correlator: 



Koz (x,x') 



K 



2 exp 



|x-x'| 



47r62 



X — X' 



The corresponding probability distribution reads: 



dWoz V (x) , dU (x)] = exp (^-^ y (x) + \? (Vf/ (x))'] dx 



(6) 



(7) 



n 



1/2 



rff/(q) 
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where U (q) is the Fourier-transformation of U (x) and Af2q is the element of volume in 
the momentum space. At 6 = Eq. ([7]) turns into Eq. ([5]). The OZ distribution has a 
characteristic energy scale Uq = k/It'/'^ . 

The disorder correlator may be regular at coincident coordinates x = x'. Then it has the 
following form: 

i^(x,x') = (f/^>/.(^^), (8) 

where (f/^) = {U^ (x)) is the average quadratic fluctuation of the random field at a point. 
The function h {u) which determines the shape of the correlation function is normalized 
by the condition /i(0) = 1. In this case the Fourier- component of the inverse correlator is 
a growing function of the wave-vector. As a consequence, (x, x') is not well defined. 
Instead we can use the probability in the momentum space: 

,W [iJ (.) , e.p (- / ^^f^..) (0) 

Simple examples of non-singular correlators are the Lorenz distribution with Hl {u) = 
(n^ + 1) ^ and Kl{^) = 27r^ (f/^) exp (— 6g) and the Gaussian distribution with 
hciu) = exp(-MV2) and i^c (q) = (271 f {U'^) exp {-b'^q'^ / 2). In all these cases the 
scale of energy is established by Uq = a/ (f/^) and the scale of length is b. 



III. SINGLE-PARTICLE LEVELS IN AN UNCORRELATED RANDOM POTEN- 
TIAL 



In this section we study the properties of single-particle states in the uncorrelated random 



potential 
Lifshitz 



defined by Eqs. (jH E])- This is a modest extension of the well-known works by I.M. 
Zittartz and Langer j^O] and Halperin and Lax 41|. Some simplifications allow 



us to find more detailed information about the distribution functions of the energy, sizes, 
distances between the states and the tunneling amplitudes between them. We are interested 
in a statistical description of the spectrum and wave functions of Schrodinger equation in 
the random potential U (x): 



^VV + (i?-f/(x))^ 
2m 



(10) 
Its energy levels E [U (x)] in a finite volume are functionals of the potential U (x). The 
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only characteristic of the random potential k together with world constants h and m estab- 
lishes the scale of length: 

(11) 



which will be called Larkin length in analogy with the scale found in the Larkin's work [l9|] 
for an elastic medium in a random field. We will show that C sets the scale of the extension 
of the deeply localized states with large by modulus negative energy. For delocalized states 
with large positive energy the same value with precision of numerical factor is the mean free 
path. The scale of energy £ can be obtained from the Larkin length in an obvious way: 



S 4 

m K 



(12) 



Further we will work in a rather rough approximation similar to that used by Larkin and 
Ovchinnikov 52] and Imry and Ma 531]. However, we start with a rigorous statement of 
the problem which gives a clue for our further estimates. The most easily calculable value 
is the density of state u {E, Q) which can be written as a path integral (see the cited works 



39|,|40|,|4l|): 



u{E,Q) = - J 6{E-E[U{^)])dW[U{^),Q], (13) 

where E [U (x)] is the spectrum of eigenvalues of the Schrodinger equation ffTOj) in the volume 
Q. 

In a large 3d volume the states with energy E ^ £ are delocalized, whereas the states 
with negative energy sufficiently large by modulus E < and \E\ ^ £ are strongly localized. 
The threshold of localization is a positive energy of the order of £. Note that this is correct 



only for 3d systems. As it was conjectured in the work by Abrahams et al. [5J], all single- 
particle states in 1 and 2d systems are localized. Near the threshold the wave function has 
very complicated fractal structure 55|], 56|]. In the interval between £^ and the transition 
from fractal to strongly localized states evolves. The latter are supported by rare fluctuations 
of the random potential, which form a potential well sufficiently deep to have the negative 
energy E as its only bound state. The problem of the deep random levels was considered 

" " n 



by already mentioned authors 



39|], [40|, [41| about 40 years ago. We reproduce some of 



their results and extend them to find the probability distribution of the levels with energy 
less than some E {E < 0; \E\ ^ £), the distances between such states and the tunneling 
amplitude between them. 
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As it is clearly seen from equations (fT][T3i) . the main exponential factor in the density 
of state can be found by minimization of the integral U"^ (x) dVL at a fixed value of the 
energy level E \U (x)], which is a functional of the random potential U (x). The latter can 
be determined as a minimum of the energy over the wave function: 



E [U (x)] = min E [U (x) , ^ (x)] = min 
Thus, we need to minimize a functional: 



F [U (x) , (x)] = J u^x)dn-xj 



t/(x) i^(x)r 



(14) 



r ft2 



|^|V^(x)|^ + f/(x)|^(x)|^ 



dQ 



(15) 



over ip (x) and U (x). Here A is a Lagrangian factor. The minimization over (x) leads to 
Schrodinger equation ( |TOl) . whereas the minimization over U (x) results in a relation between 
U (x) and ip (x): 

f/(x) = A|V'(x)|2 (16) 

Thus, equation ffTOl) turns into the Ginzburg-Landau equation. For our purpose the most 
important consequence of the relationship ( fT6l) is that the fluctuation potential well U (x) 
has the same characteristic linear size R as the wave function ipipi). It is clear that the 
maximum probability requires the bound state with the fixed energy E to be the only 
bound state in the potential well. Otherwise, at the same energy, we need a deeper well, i.e. 
larger f/^ (x). For the same reason the fluctuation well must have the spherical shape. Let 
the radius of the well is R. Then the depth of the well can be estimated as f/min 
the energy level in it differs by a factor of the order of 1/2: ~ ^2rnW^ '^^^ exponential 
factor in the density of state reads: 



exp 



2fi:2 



exp 



R 



exp 



1^1 



1/2" 



(17) 



where we redefined the Larkin length C incorporating the factor 27r/3 in it. Here and 
henceforth we perform calculations for 3d systems. The results for other dimensions will 
be derived later. The result ( fT7|) is valid provided the number in the exponent is large, i.e. 
i? < £ and \E\ > S. 

Let us consider the probability or the part of volume q (R) occupied by the wells with 
the radius less than R or the energy less than E = — 2mf>2 ■ It is obvious that q (R) contains 
an exponential factor (fT7|) and some preexponent. Since q (R) is a dimensionless value, the 
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preexponent must be a function of the dimensionless ratio 4: 



?(^) = /(4Vxpf-4) (18) 



RJ \ R 

The function / [x) must be much slower function of its argument than the exponent. Most 
naturally it is a power function / (x) ~ with a critical exponent a. It can be extracted 
from the work by J. Cardy 42[: a = 1 for the uncorrelated disorder. We will see that it 
is inessential for further conclusions. By knowledge of q (R) we can calculate the number 
nm{R) of the wells with the radius less than R per unit volume, i.e. the density of such wells 



58( 1 . In order to do that the unit volume must be divided into the cells of the volume R^. 



Each cell can contain or do not contain the fluctuation, but fluctuations with the centers 
approaching each to other at a distance less than R must be considered as one asymmetric 
potential well. The number of such cells in unit volume is R~^. Thus: 

(R) = R-\ (R) = R~'f (^^^ exp (^-^^ (19) 

The average distance d {R) between the wells with the radius less than R reads: 

d (R) = {n^ {R)r'/' = Rfl^ exp [j^ (20) 

This equation shows that the distance between the wells is significantly larger than the size 
of the wells. The situation is depicted in Fig. [21 Finally we find the tunneling factor t (R) for 
two typical wells with the energy levels of the same order or the radius R of the same order 
of magnitude. It is given by the semi classical expression t (R) = exp (— ^ J |p| c//) , where the 



integral path connects two wells. By the order of magnitude \p\ ~ ^J'lm \E\ h,/ R and the 
length of the path of integration is ~ ci [R). Thus, ^ j \p\dl d/ R = /^/^ exp (^) . Finally 
we find: 



t (R) = exp 



-/^/'exp(-^ 



(21) 



3R^ 

At i? ~ £/3 or £■ ~ —9S, the distances between the optimal potential wells become of the 
same order of magnitude as their size R. Simultaneously the tunneling amplitude between 
the wells becomes of the order of 1. The potential wells percolate and tunneling is not small, 
but the states still are not propagating due to the Anderson localization 30 1. 
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Figure 2: Deeply localized states. For low densities of bosons the distance between deeply localized 
states of radius R <^ C {or smaller) are separated by distances d{R) » C. 



IV. BOSE GAS IN A LARGE BOX WITH AN UNCORRELATED RANDOM PO- 
TENTIAL 

In the ground state of an ideal Bose gas in a large box with the Gaussian random potential 
all particles are located at the deepest fluctuation level. In the box of the cubic shape with 
the side L the deepest level which occurs with probability of the order of 1 has the radius 
R determined by equation: L^riw (R) = 1, i.e. R ~ sin{L/c) • ^^^^ that the prefactor 
introduces a correction to the denominator of the order of In (in ^) which can be neglected. 
The corresponding energy is ~ —9S (in^)^. Such a state is highly non-ergodic since 
the location and the depth of the deepest level strongly depends on a speciflc realization 
of the disordered potential. Therefore the average energy per particle and other properties 
averaged over the ensemble has nothing in common with the properties of a speciflc sample. 
Even an infinitely small interaction immediately changes the situation: the system becomes 
ergodic in the thermodynamic limit, i.e. when first the size of the system grows to infinity 
and then the interaction goes to zero. For example, the energy per particle in sufficiently 
large volume coincides with its average over the ensemble. The reason of such a sharp change 
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is that, at any small but finite interaction, the particles can not more fill one well since their 
repulsion finally overcomes the attraction to the potential well. They will be redistributed 
over multiple wells. Since the distribution of wells in different parts of sufficiently large 
volume passes all possible random configurations with proper ensemble probabilities, the 
ergodicity is established. In this section we find how the interacting particles eventually 
fills deep single-particle levels simultaneously changing their size and shape by their self- 



consistent field. Similar idea was considered in the work by Lugan et a/. [32|. However, in 
this work the consideration was done at fixed chemical potential instead of fixed number 
of particles. Besides, they did not perform geometrical analysis and calculation of the 
tunneling amplitude for single-particle states. Therefore, their results do not contain our 
explicit expressions for the geometrical structure and tunneling of the localized states as 
well as the analysis of the transition line. 

In a real experiment the Bose gas may be quenched in a metastable state depending on 
the cooling rate and other non-thermodynamic factors. This is what M.P.A. Fisher et al. 
jiol l call the Bose glass. Such a state is also possible in the case of weakly repulsive Bose gas. 
It looks plausible that the Bose glass whose metastable state has energy close to that of the 
ground state has also close to the ground state statistical characteristics, but this question 
needs more thorough study. However, as it will be demonstrated later, in the case of cooled 
alcali atoms the tunneling amplitude still remains large enough to ensure the relaxation to 
the equilibrium state in 10~^^10~^s. Our further estimates relate to the real ground state. 

The weakly repulsive Bose gas in a random potential is described by the well-known 
Hamiltonian: 

^ = E l^^i^P + I / (^^^)' / ^ (22) 

p J J 

where il) (x) = Qr^^"^ Sp'^p^^P (^P^/^) is the secondary quantized wave function; the posi- 
tive coupling constant g is associated with the scattering amplitude a by the relationship 57l |: 

(23) 



m 



As in the Bogoliubov's theory 59|] and its extension by Belyaev 60|], we assume that the gas 
criterion na^ <C 1 is satisfied. Here n = N/Q is the average particle density and is their 
total number. The purpose of the following consideration is to find how, with the growth 
of the number of particles, they eventually fill the fiuctuation potential wells. Implicitly 
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our considerations takes in account the change of the optimal potential wells due to the 
interaction. 

Let particles with the average density of particles n fill all potential wells with the 
radii less than R in the ground state. The average number of particles per well is 

'''' 

The local density inside the well of the linear size of R is np{R) = The gain of 

energy per particle due to the random potential is E (R) = — ^^^2 ] the repulsion energy due 
to interaction is equal to gnp{R) = where we used Eq. ( 123|) and the well-known 



relation for an effective potential field induced by a gas of scatterers STj. Minimizing the 
total energy per particle Etot (R) = ~ 2mR^ + j^f ^" over R and employing Eqs. (IT^ and 
(12^ . we find the value of R corresponding to the minimum of energy at fixed n with the 
logarithmic precision: 



C 
1^ 

where the critical density ric is defined as follows: 



^W = Tni7, (25) 



n 



{3C^ay' (26) 



The factor / in Eq. (fT9|) as well as the next approximation to the solution (125|) leads to the 
corrections of the type In (in ^) which can be neglected. The distances between the filled 
wells according to Eq. (l20l) are: 

din) = A:(^y^' (27) 



In ^ V n 

n 

This distance strongly exceeds the average size of the potential well (l25l) at n <^ Uc- At the 
same condition the chemical potential of atoms can be estimated as: 

a (n) = — — = In — 28) 

^ ^ ^ 2mi?2 (n) 2 V n / ^ ^ 

The tunneling amplitude t (n) between two wells separated by a typical distance r (n) can 
be found by employing the single particle result (12T|) : 



t (n) = exp 



n 



(29) 



Thus, the Bose gas at n -C n,, is fragmented into a multiple clusters of small size R (n) 
separated by much larger distances d{n) and containing about C/ 3a(ln^)^ particles 
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each. The amphtude of tunnehng between the wells is exponentially small, so that the num- 
ber of particles in each cluster is well defined. The phase is therefore completely uncertain. 
Such a state is a singlet with inhomogeneously distributed particles, a random singlet: the 
ground state is non-degenerate. 

The previous approach is valid if A/'(-R) ^ 1. From Eqs. (IT^ [2H ESD we find A/'(-R) = 
Lj 3a (in^)^ . Thus, the disorder must be weak enough to satisfy the constraint L ^3a. 
But it is not sufficient. Additionally the average density n must be large enough. Namely, 

1 /3~ 

it must satisfy strong inequalities exp — (;^) ric -C n <^ ric. At smaller n, only few 

1 /3~ 

particles appear in a potential well. At n <^ exp — (^) ric, the size of filled potential 
wells approaches 36a. Most of the filled wells contain 2 particles, but only a small part of the 
total number of wells of this size (per unit volume) (36a) is filled, approximately 0.5 x 
10~^a~^exp (— 3£j-)- The consistency of our approximation requires that the gas parameter 
remains small inside an optimal potential well. The density inside the well is equal to n^. 
Therefore, the interaction is weak inside the well if n^a^ = A <^ 1. 



3£ 



The compressibility = ^ In is finite as expected for the Bose glass phase 16 1. 

In the absence of interaction the probability that two single-particle eigenstates localized in 
different potential wells have the same energy is zero. Thus the zero temperature conductiv- 
ity vanishes in the thermodynamic limit. If the interaction is switched on, we can consider 
its effect as splitting of the energy levels by an amount of the order gn^. For n ^ iic 
the wave functions are still localized and the conductivity at zero temperature vanishes. 
Thus, the Bose glass phase is insulating. At finite temperatures however the energy gap 
between different localized states can be bridged by thermal activation leading to the Mott 



variable range hopping (VRH) 51|- Below we apply the VRH argument to our situation. 
The tunneling probability between two wells, each of extension i?, separated by the distance 
L, is I t{L,R) 1= exp(— 2L/_R) (see analogous derivation in section 3). Together with the 
activation factor it gives the hopping probability between the two wells equal to 

P(T) ~ e-2^/^e-^^/^, (30) 

where AE is the energy difference between the two distant wells. There is with high 
probability at least one localized state in the interval of energy AE provided the condi- 
tion AEv{E)E^ ^ 1 is satisfied. Employing the derived dependencies E —fi^/{2mB?), 
R = C/ \n{nc/n) and i^{E) = £~^R~^ exp {—C/R) and maximizing the probability (|30|) over 
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L, we find the hopping conductivity ct{T) 

a{T) ~ e-^[^"^/(^")l''', (31) 

where C is a constant or slow function of \\inc/n. In a similar way it is possible to calculate 
the non-linear field-dependent conductivity. 



V. TRANSITION REGION 



Though our results are valid only for n <^ rZc, it follows from them that at n ~ ric the 
overlapping of different occupied wells becomes large and the tunneling amplitudes reach 
the value of the order of 1. Therefore, the phase coherence in different wells grows until 
at some critical value of density which we roughly identify with rtc the global coherence 
is established. Thus, aX n Uc, the quantum phase transition from localized singlet to 
the superfluid state proceeds. As usual, the order parameter is the average phase factor 
{exp{iip)). The superfluid state bears clear traces of disorder: its superfluid density is 
inhomogeneous at characteristic scale C. It can be treated as a disordered superfluid. An 
alternative argument goes as follows. The condition n ~ can be rewritten in the form 
C ^ ^ where ^ = (an)"^/^ denotes the superfluid healing length. The latter describes the 
range over which the superfluid order parameter, when changed by a local perturbation (e.g. 
by a wall or in the center of a vortex), regains its bulk value. Since C denotes the scale at 
which the disorder becomes relevant, the condition ^ ^ C simply means that the superfluid 
order parameter cannot adapt to the rapid variation of the disorder which happens at the 
scale £ and superfluidity is destroyed. In the opposite case n ^ Uc, i.e. ^ ^ £ this is not 
longer the case and superfluidity is restored. 

At n ^ Uc the energy of the repulsion gn becomes much larger than the characteristic 
disorder energy S. Then the gas becomes almost homogeneous with precision of small 
parameter S/gn ^ Uc/n. In this case the approach by Huang and Meng 21[ based on the 



homogeneous zero-approximation ground state and refined by Giorgini et al. [22] and by 
Lopatin and Vinokur 23| is justified. They demonstrated that the superfluid density is 
not equal to the total density as it is without disorder, but remains close to the total density. 
Thus, the superfluidity definitely strives in the limit of large density. We conjecture that 
Us vanishes at the same quantum phase transition point n = ric aX which the coherence 
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disappears. 

The small corrections found in the cited work 



2l| can be estimated by the order of 



magnitude without long calculations. Indeed, the disorder Hamiltonian reads: 

Hdis = ^ ^qflpflp+q, (32) 

p>q 

where is Fourier transform of the random field U (x). It produces the change of energy 
in the second order of the perturbation theory. The change of energy per particle can be 
estimated as follows: 



A£ = / '-^^ (33) 



K 



.2 



The delta-like correlation function of the random potential (jl]) implies that (U^) 
for any q. The excitation energy 5q can be estimated as in the effective range of 

integration q < ^mgn after subtraction of the contribution of disorder to the energy of the 
normal state. Collecting all these factors, we find: 

With precision of a numerical coefficient 0.85 this result coincides with that obtained by 



Huang and Meng (Equation (7) of their work 



gives the value of energy per unit volume; 



it must be divided by n to find the energy of the ground state per particle). In terms of the 
Larkin length and energy and critical density introduced earlier (see Eqs. ( ITT| [T2l 1261) ) the 
correction to the energy of the ground state per particle can be expressed as 

The change of superfluid and condensate densities can be found from their initial values 
(one exactly and another approximately equal to the total density n) by multiplication to a 
small factor proportional to Ae/eo ^ ^/C, where Eq = ^ = ~ 

In discussing critical phenomena it is crucial to determine the critical (marginal) dimen- 
sions. In order to do that we need some results which will be derived later. We will show 
that the critical density at the transition between the superfluid and the Bose glass phase 
in d dimensions reads 



'^^ ~ nd-2 /-2 ' ~ ^d, (36) 
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where Cd = (^^/m/t)^/^^"'^) and ~ (a'^~^?7,)~^/^ denote the d-dimensional Larkin and 
heahng length, respectively. For strong disorder and strong interaction the natural limits of 
the validity of equation (136!) seem to be ^ n~^^'^ and Cd ^ n~^^^, i.e. the healing length 
and the Larkin length must be much larger than the mean spacing between the bosons. 

For weak disorder, the Larkin length diverges at — ^ 4. Hence in dimensions d > 4 
disorder has to overcome a threshold value to destroy superfluidity even in the limit of 
weak interaction. Thus c? = 4 is an upper critical dimension of the problem as has been 
discussed already in A lower critical dimension is = L This can be seen from a 

calculation of fl35|) in one dimension which is logarithmically infrared divergent reflecting 
strong quantum fluctuation. Giamarchi and Schulz [l8| in a seminal paper found a quantum 
phase transition in Id Bose gas in disordered environment from a superfluid to an insulating 
phase for increasing interaction parameter K ^ 2/3 (compare Appendix). Since then this 
transition has been considered as the genuine transition between these two phases (see e.g. 

n 

|16l|). However, at this transition ^in ^ 0.25 strongly violating the gas condition C,in ^ 1. 



Thus, we conclude that the transition found in [18|] is not the one-dimensional equivalent of 
the transition taking place at n ~ ric- The phase portrait obtained in this work is reproduced 
in Fig. [HI It follows from this figure that there are two different phase transitions at small 
disorder: one at large value of dimensionless coupling constant K and another at small K. 
The large part of the phase diagram beyond a small vicinity of the critical point K = 2/3, 
K = was obtained by a speculation rather than a rigorous treatment. Therefore, the 
question about the existence of two different quantum phase transition remains open. 

Critical behavior near th e q uantum transition was studied by Halperin, Lee and Mafisl 
and by M.P.A. Fisher et al. 16[. In the latter work the scaling consideration was developed 
and even the critical exponents were found. However, we are not sure that they are valid 
for our problem since, according to previous remark, they are associated with the strong- 
coupling critical point. Secondly, even if there is only one critical point, the marginal 
dimensionality is 1-1-1 and it is very doubtful that exact critical exponent can be found in 
the case 3-1-1. 
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VI. CORRELATED DISORDER 



In this section we consider properties of the weakly interacting Bose gas in a big box 
with a correlated random potential. The properties of correlated random potential were 
already described in the second half of Section [ITl (see the text related to Eqs. ([Hl-in]))- Here, 
following the line of consideration developed in Sections IllUlVt we study first the size and 
distribution of single-particle states and tunneling between them and how they are modified 
by interaction. We start with the single-particle states. 



A. Single-particle states in a correlated random potential 



Keldysh and Proshko 



431, and Kane 4^ were the first to find the electron density of 



states in a semiclassical random potential. They proved that, in contrast to the uncorrelated 



disorder for which v [E) ~ exp — [E / £) 



.1/2 



, the density of state in a correlated disorder is 



a Gaussian function of energy v {E) ~ exp \^—c{E /Un)^\, where c is a numerical constant, 
ndependently but later the same result was found in 45|, |46 



47l |. Shklovskii and E: 
48l | derived the same result employing the instanton approach. John and Stephen 



ros 



49| 



found the dependence of the preexponential factor on energy. Their result was confirmed 
by Thirumalai [s^ by a different method. Below we rederive their result and find new 
important characteristics of optimal fluctuation wells in the case when the correlation length 
is large enough. It is intuitively clear that, at 6 <^ £, the density of state and other 
characteristics of the spectrum and optional states only slightly differ from their values for 
uncorrelated random potential. At 6 ~ £ numerical coefficients will be different from those 
for uncorrelated disorder, but with this reservation still our semiquantitative description is 
valid. Therefore, the most interesting is the situation with b ^ C Such a disorder we call 
strongly correlated. In what follows we assume the strongly correlated disorder . The Larkin 
length for the correlated disorder reads: 



(37) 



It is convenient to introduce another characteristic length B = (^)^^^ (nSo 
strongly correlated disorder B <^ h. The Larkin length is associated with B by the fol- 



1/2 



For 
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U{x) 




Figure 3: Long range correlated potential. There are typically many bound states per well. 



(38) 



lowing relation: 

C _ /B 

To find the exponent in the density of deeply localized states it is necessary to minimize 
the functional J U (x) (x, x') U (x') rfxrfx' at fixed energy of the quantum state in the 
potential U (x). The condition of minimum is reduced to a following relationship between 
the optimal fluctuation potential U (x) and the wave function ip (x): 

f/ (x) = A j K (x, x') (x') 1^ rfx', (39) 

analogues to Eq. (fT6l) for the uncorrelated disorder. This equation shows that the char- 
acteristic size of U (x) is always equal to h. Indeed, it is correct if the characteristic size 
of the wave function less or equal to h. In the opposite case the characteristic size of the 
potential would be much larger than h and, respectively, the probability of such configura- 
tion would be much smaller. Deep levels have negative energy much larger by modulus 
than the energy scale of the random potential Uq. Since the optimal fluctuation potential 
well must have the level its depth is not less than l^'l. The fundamental difference 
between the optimal potential wells in the cases of uncorrelated and strongly correlated 
disorder is that the former contains only one level, whereas the latter contains many lev- 
els as shown in Fig [31 Indeed, the number of levels in the optimal well is not less than 
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^ (^^2m\E\h/h)^ ^ ^ {^fhnU'Qh/hY ~ {^f^^ = (|)^ > 1. Therefore, the level E is 
located close to the bottom of the optimal potential well. Its depth with high precision 
is equal to \E\ and its shape is determined by the function h{r/b) from Eq. (JHl). If the 
correlation function is non-singular, then U (x.) = E h{r/h). The exponent in Eq. ([T]) can 
be calculated exactly leading to the result for the probability q (E) to find the level with 
energy less than E (with precision of a prefactor): 

q{E)=exp(--^) (40) 



2(f/2) 

In a more accurate treatment one needs to take in account that the level E can be not the 
ground state in the optimal well, but an excited state. Since even excited states are close to 
the bottom, the spectrum is the same as for a spherically symmetric quantum oscillator with 

1 /2 

the frequency lo = {\E\ h" (0) /mlP') . Therefore, the depth V of the well is not exactly 
equal to \E\. It depends on what excited state occupies the particle: 

Ki,n2,n3 = 1^1 + [(^1 +n2 + ns) + 3/2] nw (41) 

The second term in the depth of well determined by this equation is much less than the 
first one. After summation over all rij (2 = 1,2, 3) from to infinity we find the corrected 
expression for the probability q (E): 



q (E) = exp 



^2 + 31^1^ 



exp 



\E\ Hio 



(42) 



2([/2) 

The ratio \E\ hwj (f/^) which appears in Eq. (142!) by the order of magnitude is equal to the 
product of two ratios i\E\ jU^f^'^ x {B /h). The first factor is large, the second one is small 
and their product varies from small to large values always remaining larger than [B/h). 
Nevertheless, a strong enhancement of probability up to a factor (b/B)"^ in comparison to 
the initial expression (j^Oj) is reached at decreasing modulus of energy. Summarizing, we 
conclude that the participation of excited states may change a comparatively slow-varying 
prefactor leaving the principal exponent ffTSjl unchanged. 

Next we analyze the case of a singular correlator using the example of the Ornstein-Zernike 
distribution. The relationship fl5I?]) is valid for singular correlators implying that the radial 
size of the optimal potential well is equal to the correlation length b. We use the functional 
([7]) to estimate the exponent in the probability q{E). The two terms in the exponent ([7]) 
proportional to U"^ and b^ (Vt/)^ are equal. With the precision of a numerical constant 
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-r/6 



the potential is equal to E {h/r) e 
accordance with the earlier results 49|], 



he integration leads to q{E) = exp {—cE^/Uq) in 



501]. Further we incorporate the constant c into the 



definition of Uq and write the equation for the probability in all cases as: 

g(i?)=exp(-EV2[/2) 



(43) 



For simplicity we put the prefactor to be equal to 1. 

The density of potential wells containing deep levels lower than E for strongly correlated 
random potential reads: 

{E, h) = exp {-E^l2Ul) (44) 
and the average distance between such wells is 

d {E, b) = b exp (EyeU^) (45) 

The tunneling amplitude between such wells can be found in semiclassical approximation: 

^/2m\E\d 



t {E, b) = exp 



(46) 



exp 



--^exp (EV6f/o^) 



It is always very small, but the direct tunneling between wells does not play an important 

ity edge in contrast to 



351 ] ■ the percolation of 



role since the particles remain almost classical even near the mobi 
the uncorrelated disorder. According to the Shklovsky estimate 
wells happens at E = — 0.9f/o. The characteristic De-Broglie wavelength H/ \phn\E\ at 
this energy is much less than b for strongly correlated disorder. Therefore, the particles 
propagate according to classical mechanics. The mobility edge energy is negative and not 
far from the percolation threshold. 



B. Weakly repulsive Bose gas in strongly correlated disordered potential 

In this section we consider the Bose gas with the average density n. The gas partly fills 
all potential wells with the depth less than E. Let the particles fill a typical well up to the 
radius R < b. If the correlator is regular at x = x', the following equation gives a rough 
estimate of energy per particle: 

, (E, R-^b,n)^E^l--y^ (- j exp j (47) 
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The first term in this expression is the gain of energy due to the random potential, the 
second term is the energy of interaction. The quadratic term in the brackets corresponds 
to the oscillator-like potential at small R. Minimizing n {E, b, R) over E and the ratio R/b, 
one arrives at following results: 



E^-Uo^2ln'^- (48) 
The critical density ric is defined as follows: 

All results are valid at n <^ Uc- According to Eq. poj) . the ratio is small. This fact justifies 
the approximation of quadratic potential in Eq. ( H71) . It allows to estimate the characteristic 
momentum of particles or the width of their distribution over momentum: 



n 



-1/2 



Ap = mujR = mujb ( In — 1 (51) 



The chemical potential is readily derived from Eqs. 



/i(&,n) = -f/oj21n^ 

n 



1-^fln!^ 
4 V n 



Two filled wells are separated by a typical distance d {b, n) much larger than the size of the 
well: 

d{b,n) = b(—y^^ (52) 



(53) 



n 

With this value of d it is ready to estimate the tunneling coefficient: 

b /n,\i/3l 



t (b.n) ~ exp 



B \ n 



The number of particles in each typical well can be estimated as J\f = ^ = nj)^ ^ 
It is large if 6 > a. 

These results conclude the description of localized states for strongly correlated disorder. 
At energies between —Uq and Uq the motion of particle is strongly disordered, but it becomes 
more and more free when energy approaches and exceeds Uq. At E ^ Uq, the random 
potential can be considered as a perturbation. The quantum phase transition from the 
localized random singlet state to a superfiuid proceeds at = Uc- Note that the critical 
density in the presence of strongly correlated disorder does not depend on the correlation 
length b. 
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VII. BOSE GAS IN A TRAP WITH DISORDERED POTENTIAL 



In this section we consider deeply localized states in a harmonic trap supplied with a 
Gaussian random potential. The Hamiltonian of the system has an additional term: the 
potential energy of Bose particles in the trap: 

Htrap = j Vtrap (x) (x) ^ (x) C?X, (54) 

where the harmonic potential of the trap generally has a form: 

Vtrap (X) = y + + U^y) (55) 

In this section we consider only the isotropic trap with Strongly 
anisotropic traps can be considered as a system with reduced dimension 2 or 1 and will 
be considered in the next section. A new scale of length associated with the trap is the 
well-known oscillator length I = \/fi/ (muj). 

A principal difference between the Bose gas in a box and in a trap is that in the latter case 
the gas forms a cloud whose size is determined by energy minimization at a fixed number of 
particles, whereas in the former case the size of the box and, therefore, its average density 
n is fixed. Nevertheless, we will see that all phase states of the Bose gas in the box appear 
when this gas is placed into a harmonic trap, but their appearance is regulated by the total 
number of particles in the gas N instead of their average density. 



A. Trap with uncorrelated disorder 

Four competing parts of a Bose-particle energy are: kinetic energy, the confining potential 
energy of the trap, the repulsion from other particles and the energy of the random potential. 
Two of them, the interaction with the trap and the random potential tend to confine and 
localize the particle. Let the gas form a single cloud of the radius R. The rough estimate of 



the two confining energy contributions can be done as follows 



34| 



\^trap - - — ^ (56) 



The disorder energy for the same cloud can be estimated as 



y-.. = —T^'--l:(iY' (57) 
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The latter estimate can be justified as follows. The cloud chooses such a position that the 
random potential energy is negative. A characteristic fluctuation of energy in a volume Q is 



1 



U fx) rfx 



1/2 



1 
o 



K (x, x') d'x.d'x! 



1/2 



Substituting into the last expression Vt = ^nR^, we arrive at Eq. ( 1571) . Comparing the two 
confining contributions, we conclude that, at C ^ i, the influence of disorder is negligible. 
The remaining two terms of energy: kinetic energy 



(58) 



and the repulsion energy per particle: 

Vint 



gn 



?>Na 



(59) 



m 

are deconfining. The kinetic energy dominates at i? ^ 3A^a. In the opposite case R <^ ZNa 
the repulsion energy dominates. Below we analyze several limiting cases. 



1 Weak disorder: C ^ i. Depending on the relative value of the interaction one should 
distinguish two different situations: 

la. Weak interaction: 3Na <^ i. In this case the interaction can be neglected. Minimizing 
too remaining terms, the kinetic energy and energy of the trap, we find R = £. Physically 
it means that all particles are condensed at the oscillator ground state, 
lb. Strong interaction: 3Na ^ £. Neglecting kinetic energy and minimizing remaining 
energy of traps plus the repulsion energy, one finds the result known as Thomas-Fermi 
approximation 61|: R= (^^Nai'^Y^^ ■ 

More interesting is the case of strong disorder. 
2. Strong disorder: C <^ i. Again we consider several different limiting cases depending 
on relative strength of disorder, interaction and the trap potential. 

2a. Weak interaction: 3Na <^ C <^ i. In this range of variables the non-ergodic phase is 
realized. Since interaction is negligible, the particles find a random potential well with the 
deepest level and fall into it. Let such a well can be found at a distance ~ L from the trap 
center. Its depth typically is about 9£ Iv? {L/C). This gain of energy must be not less than 
the loss of the trap energy ma;^L^/2. A typical value of L appears when both this energies 
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Figure 4: Regime diagram of atoms in three dimensional traps: uncorrelated disorder. K denotes 
the size of the single existing atomic cloud. L is the size of the cloud of fragments. 



have the same order of mae; nitude. Then L ^ 6^2 {H^ I C) In (£/£). A typical size of the well 
isi?^/:/ (6 In (£/£)). 

2b. Moderate interaction: L ^ 3A^a <C I. In this case the ergodicity is restored. Our 
experience with the gas in a box prompts that the gas cloud is split into fragments each 
occupying a random potential well from very small size till some size R depending on A^. 
From Eq. ( l28l) the typical disorder energy per particle is /i = —E (in^)^. It becomes 
equal to the trap energy at a distance L ~ (^^/£)lnr (see the definition of the parameter 
F below in Eq. (l60l) ). Therefore, the average density is n ~ N C? ji^ . The ratio n^jn plays 
the decisive role: the state of the Bose gas is fragmented and strongly localized when it 
is large; the transition to delocalized superfluid state proceeds when this ratio becomes 1. 
With precision of a numerical factor the ratio n is equal to a new dimensional parameter: 

Strongly localized fragmented state is realized at F ^ 1. The transition proceeds at F = 1. 



The superfluid phase emerges at F < 1. The regime diagram is shown in Fig. HI 6J] Note 



the counter-intuitive dependence of the cloud size on the number of particles: the cloud 
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slightly contracts with increasing number of particles. It happens because the number of 
particles in each fragment increases more rapidly with the average density than the number 
of fragments. 

B. Trap with strongly correlated disorder 

Following the line of development accepted in the previous subsection, we start with the 
case of weak disorder. 

1. Weak disorder: -f^fUo fku or B ^ i. The value in the left-hand side of this 
interaction is the characteristic variation of the random potential on the minimal of two 
linear sizes: oscillator length i or correlation length b. In this case the random potential 
only slightly changes the state of the gas in the absence of disorder (see previous subsection). 
In what follows we assume for simplicity that £ > b. 

2. Strong disorder: Uq ^ tiw or B <^ i. As in the previous subsection we consider 
different relative strength of the interaction and disorder. 

2a. Weak interaction: 3Na <^ i. The interaction is negligible. Therefore, all particles 

find the deepest well in the trap and fall into it. The typical distance from the center of the 
trap to the deepest well is L = ^ (\n^^ and its typical depth is E = —Uq^Q In ^. In 
this range of parameters the system is non-ergodic. 

2b. Moderate interaction: ^ <^ 3Na <^ £. In this case the particles are distributed 
among multiple potential wells. The system is ergodic. The average depth of the well is 
E (N) = -UoV2 InT, where 

7-/5/2 .6 

r = = — fell 

is the dimensionless parameter analogues to that introduced in the previous subsection. It 
represents the ratio Uc/n. Until F is large the ground state of the Bose gas is deeply localized 
random singlet. The transition to the superfluid state proceeds at F ~ 1. The total linear 
size of the cloud is L = 2^/^^ (InF)^^^. Inside a fluctuation potential which has a typical 
size b particles occupy a smaller sphere of the radius R ~ [Nab'^B'^Y^^ . The phase state 
(regime) of the system depends on three dimensionless parameters: B/b, 3Na/i and i/b. 
The regime diagram is shown on Fig. [5l 
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Figure 5: Regime diagram of atoms in three dimensional traps: correlated disorder. R denotes the 
size of the single atomic cloud, L denotes the size of the fragmented state. 

VIII. LOWER DIMENSIONS 



In the experiments with cooled gases the trap is mostly realized as a rotation ellipsoid 



with a large aspect ratio. Therefore, and also for comparison with exact Id solution 
is reasonable to extend our theory for lower dimensions d = 1,2. 
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A. Uncorrelated disorder 



We start with the extension of Eq. (fTTj) connecting the Larkin length Cd and the constant 



k"^ for the uncorrelated disorder: 



(62) 



following from the dimension consideration. The characteristic energy reads: 

1 fm'^K^^ ^ 



2mCl 2 V 



(63) 
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The extension of the exponential law (fT71) for the probability to find the deep potential well 
of the size not larger than R is: 

4- 



q{R) = f exp 



R 



/exp 



\E\ 



4-d 
2 



(64) 



where the Larkin length differs from that given by definition fl62l) by the factor {0.^/2)^^ , in 
which Vtd is the volume of the sphere with the radius 1 in (i— dimensional space. As earlier 
/ is a power-like function of the ratio C/R at C/R^ 1: /(x) oc with a = 1 for d = 2,3 



and a = for d 



42]. The extension of the Eq. f|T9l) for the density of the fluctuation 



potential wells with the radius not exceeding R looks in d dimensions as follows: 

4-d" 



{R)=R-''fexp 



and the average distance between such wells is 



d{R) =Rr^/'^exp 



R 



d \ R 



(65) 



(66) 



The tunneling amplitude between nearest wells has a characteristic value 



t (i?) = exp 



-l/d 



exp 



IfCd 
d \ R 



(67) 



Next we consider the filling of the potential wells by the diluted Bose-gas of the density 
n = N/L'^. Again, let the gas fills all the wells up to the radius R. The average number of 
particles per well is: 

4-rf1 



AfiR) 



n 



riyj (R) 

and the density of particles in the well is 



/"^rfexp 



Cd 
R 



(68) 



n, 



(R) = (R) / {QdR'') 



n 



exp 



R 



(69) 



The repulsion energy per particle is as before gdrip^R), but the connection between the 
coupling constant g and the scattering amplitude a depends on dimensionality of the space: 



9d 



m 



(70) 
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At c? = 2 this equation is invalid, but the couphng constant g2 is well defined and has 
dimensionality energy x length^. Minimizing the energy per particle in complete analogy 
with what was done in Section HVl we find the maximal radius of the filled well as a function 
of density: 

R{n) = —^, (71) 
(1^ 

and the critical density is defined by the equation: 

^d^d ^d 



rir 



(72) 



In 1-d Bose-gas the critical density increases with the repulsion strength in contrast to 2- 
and 3-dimensional cases. The reason is that the 1-d Bose-gas with strong repulsion (hard 
core) is equivalent to the Fermi-gas whose kinetic energy increases with the size of the core. 
The ratio Uc/n can be rewritten in equivalent form: 

^=(1^)', (73) 

where = H/ y^2mgn is the coherence or healing length of the Bose-gas without disorder. 
The condition n -C at which the deeply localized ground state is realized implies also Cd <^ 
C,d as it could be expected. In Id case large values of the gas parameter nai correspond to the 
almost ideal case in which kinetic energy is much larger than the interaction. Accordingly, 
the gas parameter inside the fiuctuation well ricCi'^ ~ (a/Cd)'^ must be large in Id case. 

We do not show other calculations analogous to those produced in Section [IV] and present 
only the final results: 
Number of particles in a filled well: 

in) = (in^Y^' ( (74) 



Distance between the filled wells: 



n 



l/d 



d(-) = ^4%) (In^) (75) 



Chemical potential: 
Tunneling amplitude between wells: 



/i (n) = -Sd ( In ^) (76) 



t (n) = exp 



n 



n 



(77) 
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Considering the gas in a trap which has a shape of a disc or a cigar, one should express 
the parameters of effective d— dimensional problem in terms of parameters of the 3d gas: the 
particle density n, 3d-scattering length a, Larkin length £, the transverse oscillator length 
£_L and the longitudinal oscillator length £. To simplify the results we consider only the most 
interesting and experimentally accessible situation of a weak confinement: 

a < £x < (78) 
Then it is straightforward to see that the effective d— dimensional Larkin length is related 



to the 3d Larkin length £ by a following relationship 34l |: 

Cd-{Ce^^)^' (79) 

Indeed, the argument in exponent of the Gaussian distribution for uncorrelated disorder is 
~ U'^VI/2k'^. The volume Vt can be approximately factorized: Vl ^ VldR'^Vt-i^d^^^'^ ■ Using the 
standard estimate for U = we arrive at Eq. (IMl) with Cd defined by Eq. ( 1791) with 

precision of a numerical factor. To find the effective scattering length a^, one should to use 
an obvious identity: ng = UdQd (both expressions are the interaction energy per particle) 
with Ud = n£^j^'^Q,-^_d being the density in the reduced dimensionality space. From these two 
equations it follows that: 

9.= '^ (80) 
Employing Eqs. (170]) and (I23p . we arrive at following relationship: 

ad~(a^l-')^ (81) 



This equation is invalid for d = 2 and instead Eq. ( IHOl) must be used. Below we consider 
possible regimes similarly to what we did for 3d system. 

Weak disorder: Cd ^ or equivalently C ^ The random potential is negligible. 

The competition of the kinetic energy, trap energy and interaction provides two different 
regimes: 

Weak interaction: gdN <^ ft^i'^^'^/m or equivalently Na <^ In this case the 

longitudinal size of the cloud R coincides with i. 

Strong interaction: Na ^ This is the region of the Thomas-Fermi regime: R ~ 

(iVa^V^i"'^)'/^'^'^ 

Strong disorder: ^ ^ or equivalently C -C 
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Figure 6: Regime diagram of atoms in a one dimensional trap: uncorrelated disorder. In reduced 
dimensions R denotes the longitudinal size of the single atomic cloud, L denotes the longitudinal 
size of the fragmented state. 



Weak interaction: ^^^i^d ^ or equivalently Na <^ ll.^~<^i^ . This is non-ergodic situ- 
ation: all particles fall into the deepest potential well of the typical size R ~ Ldl ^In ^ j ^ . 
For (i = 2 it happens at -C — ; for (i = 1 the non-ergodic regime exists if ^ [a^ Cf^'^ . 

2 (3-d) 

Moderate interaction: L'^^t^ ~ <^ Na <^ i . In this range of interaction the 

cloud consists of multiple fragments. They are weakly coupled each to other provided the 
dimensionless parameter F = t^'^ j {Na'^f'^Cf''^) = P'^ / I NaC^£_^ is large (note 

that the exponent at i± is zero at (i = 1, 3 and -1 at (i = 2). We remind that F is equivalent to 
the ratio Uc/n. A typical longitudinal size of the fragment is i? = Cd/ (In F)^''*-^"'^-*; a typical 
distance between fragments is d = RT^^'^; the tunneling amplitude between fragments is 
t ~ exp(— F^/'^). The transition from localized state to the phase-correlated superfiuid 
state proceeds at F 1. At F ^ 1 the disorder is negligible and the Thomas- Fermi 
approximations works. The phase diagram in one dimension is shown in Fig. [6] 
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B. Strongly correlated disorder 



As in 3d case, the distribution of random potential is characterized by the correlation 

1 

length b and the characteristic energy Uq. It is assumed that h -> C [K^ / {w?U^h'^)] 
or equivalently b B = :^7=^- The characteristic size of a deep potential well is equal to 
the correlation length b. The probability to find the well containing the level with energy 
E < and sufficiently deep, so that \E\ ^ Uq, is given by the Gaussian exponent Eq. PUI) . 
The number of such wells per unit volume is 

n^,{E)=b-''exp(^-^^ (82) 

and the distance between them dd{E) = [n^diE)]^^^'^- The gas with the average density 
rid fills the wells up to the depth \E\ = Uo^yin (ricd/nd). The critical density is given by 
equation: 

The latter formula is invalid at = 2, but the former one, employing the value gd^ remains 
correct. The particles in each well fills a d-dimensional sphere of the radius R = bj^ = 
I ^ =. The transition to the superfluid state proceeds at nd = Ucd- 
Proceeding to the gas in a trap, we first specify conditions which are realistic and allow 
to reduce the number of possible regimes: 

a<^i±<^ C<^b (84) 

At these conditions fi'd^l"'^ = 9', o,d~'^^^±'^ = and Eq. (183|) can be rewritten as follows: 

_ QdUoil-" _ Qd /o^N 
""^^ - ~d^~ - A^dB^a ^^^^ 

Note, that this equation is valid at (i = 2 as well. The possible regimes are classified as 
follows. 

1. Weak disorder: Uq <^ fuu or equivalently B i. 

la. Weak interaction: Na <^ In this case the size of the cloud R = i. 

lb. Strong interaction: Na ^ The Thomas-Fermi range: 

And Nat\^ 
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Figure 7: Regime diagram of atoms in a one-dimensional trap: correlated disorder. In reduced 
dimensions R denotes the longitudinal size of the single atomic cloud, L denotes the longitudinal 
size of the fragmented state. 



2. Strong disorder: Uq ^ huj or B i. 

2a. Weak interaction: J^i^z-d ^ f^o or equivalently Na <^ . Non-ergodic situation. A 
typical depth of the deepest potential well is = — f/o\/2(iln {L/b). A typical distance of 
this well from the center of the trap is L ~ \/Uq/ (muj^) ^2(iln|^^ = ^ (^d\n^^ 

2b. Intermediate interaction: — g^— <^ Na <^ £ ij_ . This is the range of frag- 
mented state. The particles in fragments are deeply localized if the parameter F = 
£2d£3-d^ (S'^+^aA^) is large. The particles in a fragment occupy a d-dimensional sphere 
of the radius R{N) = b/V\nT. Their typical energy is E = — [/qVIu F. A typical distance 
between fragments is d (N) = bT^^^. The number of particles inside each fragment is about 
fe'^^i"'^/ {aB'^). The total size of the cloud is ~ (f/B) (InF)^/^. The transition to the su- 
perfiuid states proceeds at F ^ 1. The phase diagram in d-dimensions is shown in Fig. [7] 
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IX. DISCUSSION AND CONCLUSIONS 



The theory allows prediction of many quantitative characteristics of deeply localized 
states. It also gives simple formulas for the critical density at which the quantum phase 
transition from the random singlet to superfluid state proceeds. Below we discuss how these 
predictions can be checked in experiments with cooled atomic gases. Four important param- 
eters can be controllably and independently varied in the experiment. They are: number of 
particles A^; the frequency uo or equivalently the strength of the trap; the scattering length 
a (by approaching one of the Feshbach resonances); the strength of disorder Uq. Using this 
freedom it is feasible to pass all spectrum of regimes described above. In particular, varying 
Uq it is possible to transit from uncorrelated to strongly correlated disorder. A simple esti- 
mate shows that, at 6 ~ 1/im, the transition from uncorrelated to strongly correlated regime 
proceeds at frequency of disorder potential ujd = ^J^lXj^Jmip- ~ IkHz which is accessible. 

Simplest experiments are the measurements of the cloud size L as a function of different 
variable parameters in the regime of multiple localized fragments. Our theory predicts that 
in the regime of uncorrelated disorder the size of the cloud is inversely proportional to the 
frequency of the harmonic trap u and proportional to the square of the disorder strength. It 
also predicts a weak dependence of the size on the number of particles ~ In A^. In the case of 
strongly correlated disorder the size of the cloud is again inverse proportional to uj, but the 

1/2 

dependence on the disorder strength is much weaker L oc Uq . The dependence on A^ also 
is weaker than in the uncorrelated regime: L oc (InA^)^''^. It would be important to observe 
a crossover from non-ergodic state with one or few fragments to the ergodic state with many 
fragments and check that it happens at A^ = £/3a for uncorrelated disorder and at A^ = 
for strongly correlated disorder (see modifications of these results for = 1, 2 in Section 7). 
Counting the number of fragments, one can find the number of particles per fragment and 
compare it with theoretical prediction. Direct measurements of the size of fragments and 
distances between them are much desirable since they give an important information on the 
balance between the inter-particle repulsion and their interaction to random potential. 

Another feasible experiment is the time-of-fiight spectroscopy after switching off both 
the trap and the random potential. In this experiment the distribution of particles over 
momenta (velocities) is measured. The distribution of momentum has a finite width Ap. 
According to the uncertainty principle the value H/Ap can be treated as the average size 



35 



of the fragment R = C/lnT for the uncorrelated disorder (see the definition of F in Eq. 
f l60|) ). The fastest is the dependence of this value on the strength of disorder: Ap oc Uq. 
Others dependencies are logarithmic. For example = — ^ = 3 /j^^ • Installing a counter 
close to the trap at a distance comparable to the size of the trap would allow to register 
the oscillations of the particle flux due to discrete character of the fragmented state. This 
is another opportunity to find the distances between fragments. In the case of strongly 
correlated disorder the width of the distribution is determined by Eq. fl5T]) . in which the 
ratio ric/n must be substituted by F defined by Eq. (16T|) . In this case Ap depends linearly 
on u. Other dependencies are much weaker. Nevertheless, the prediction that the ratio 
Ap/uj depends on one dimensionless parameter F can be experimentally checked. 

The transition between localized and delocalized coherent state in the random potential 
was confirmed in several experiments (see Introduction). We propose to make more detailed 
measurement of the transition manifold and check our predictions formulated in Eqs. (l26l 

EUD. 

An important question is whether the relaxation to the ground state can be reached 
during a reasonable time interval compatible with the time of experiment. We analyze this 
question for the uncorrelated or weakly correlated disorder. In this case the relaxation time 
due to tunneling can be estimated as r = 27ruj~^t~^, where Un ^ (In F)^ is the characteristic 
frequency of the optimal potential well and t ~ exp [— F^/'^] is the tunneling coefficient (see 
Eq. ( !29l) ). For numerical estimates we accept F ~ 125, I ~ lOyum, b ^ C ^ Ifim, a ~ 
0.01/im, ^ 27,000. Then t^^ = 148 and r ~ 0.06s. The Larkin length can be increased 
by decreasing the amplitude of the random potential (£ oc Uq'^). Simultaneously at fixed 
values A^, i and a the value F decreases as C~^. This example shows that the equilibrium is 
accessible, though it may be difficult to reach large ratio C/b. In experimentally important 
case of Id (cigar-like) cloud the relaxation time is shorter since it is proportional to exp(— Fi) 
with Fi = j^^- At i = 25/im, C = l/im and a = O.Olfim the ratio Fi reaches value 1 at 
A^ ~ 85,000. 

In conclusion, we presented a simple physical picture of deeply localized states of the 
Bose gas in a random potential. We demonstrated that the particles eventually fill the deep 
potential wells formed by fiuctuations of the random potential and by their self-consistent 
field. Based on this idea the geometrical and physical properties of these states, such as the 
size of the clusters, the distances and the tunneling amplitudes between them, the size of the 
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whole cloud formed by the gas in a harmonic trap were calculated with precision of numerical 
coefficients. We discovered that the physical properties of these states depend significantly 
on correlations of the random potential. It occurs that the ground state of the system can 
be either almost homogeneous and coherent (superfiuid) if the disorder is weak enough, or 
fragmented and strongly localized. In particular, if the disorder is much stronger than the 
repulsion between particles, the system transits into a non-ergodic state, whose properties 
even in the equilibrium strongly depend on the specific sample. At growing number of 
particles the system occurs in an ergodic, but strongly localized ground state consisting 
of multiple particle clusters populating deep fiuctuation wells. At the number of particles 
increasing the tunneling between different potential wells increases leading to the phase 
correlation and finally to the quantum phase transition to the coherent (superfiuid) state. 
We have found simple expression for the gas density at this transition in the gas confined 
by a big box and equation for the phase transition manifold in the gas confined by the 
harmonic trap. It was shown that in any dimension d < A the repulsive interaction causes 
delocalization and quantum phase transition to the superfiuid state at a critical density. 
Thus, the interaction overcomes the Anderson localization of all single-particle states in 1 
and 2 dimensions. In 3 dimensions the single particle levels which are filled at critical density 
have energy still below the mobility threshold. 
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Appendix: The superfluid-insulator transition in one dimension 

In this appendix we discuss the relation of our approach to the theory of Giamarchi and 



Schulz 



It is convenient to state first the relation between the mass m, the boson density 



n, the scattering length a, the parameters K and k of the field theory of Id bosons [18|, |37 |. 
We omit all subscripts i which are excessive for one dimension. The phonon velocity v is 
then given by 

v=( 87 

37 



where x denotes the compressibilty of the bosons 

dn 1 



du TihvK 



(88) 



The Haldane inverse compressibihty >ch 37|] is related to our standard definition by x 



V? I hh- Indeed, according to our definition the elastic energy per unit length is i^^ ~ 
which implies 

1 /m\i/2 

K = — [ — ) . 89 
The Boson coupling constant K is inverse to that of fermions. These relations agree with 

n n 

those of taking in account that in pj|] the density is dimensionless and h= 1. To get the 



relation between a and K we quote Haldane [s^] referring to earlier work of Lieb 62|. Using 
the dimensionless parameter 7=1/ (ir^na) (we absorbed a factor I/tt^ in the definition of 
7), i^' is related to 7 by 

K ^ 7^/M 1 - -7^/2 ) , 7<1 (90) 



2 

^1-— j , 7»1. (91) 

with a smooth interpolation in between. Thus K changes continuously from (for a ^ 00 
corresponding to no interaction) to = 1 (for infinitely strong interaction). Since for 
K = 1 bosonic and fermionic description coincide, an infinitely strong repulsive interaction 



corresponds to free fermions 63| . Let us formulate the condition of superfiuidity in terms of 
coupling constant K. The condition of superfiuidity for Id disordered system n > Uc ^ a/C"^ 
(see Eq. f l5B]) ) can be rewritten as 

1 /mK\2/3 1 
K>K, = — ( — ) 92 
Tin \ h"^ / TvnL 

where k denotes the disorder strength. We will show that K ^ K,^ remains the phase 
boundary between the superfiuid and Bose glass phases also if quantum fiuctuations are 
taken into account. 

Let us study how quantum fiuctuations renormalize the disorder strength k. In lowest 
order in k the disorder strength is renormalized by the Debye- Waller factor: 

K,s ^ Ke-^^^'l (93) 

Here 
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Figure 8: Phase diagram in one dimension. 



where we have used the relation for small K. Thus, for K < K,^ there are no fluctuation 
corrections, and our condition fl36l) tells us that there is no superfluidity. In the opposite 
case K > K,^ we find (by making the calculation self-consistent) for max(i^'K, 2/3) < K 



"if 



(95) 



Approaching the critical value K = 2/3 the effective disorder strength is renormalized to 
zero for i^^ < K < 2/3. In the opposite case there is no renormalization and the disorder 
is relevant. From this considerations it is clear that the transition happening at K = 2/3 
occurs only in one dimension. There is another transition line K = which exists in 
higher dimensions. In the next order in /t, the condition K = 2/3 changes to = i^^ with 
K^=Q = 2/3 18|. The transition at K = i^^ is of the Kosterlitz-Thouless type. Even in a 
pure system superfluidity persists only when the healing length ^ is smaller than the system 
size L, i.e. for ^ ~ (a/nY^"^ < L. For small K this condition can be written as irnK > 1/L. 
Thus, for a oo the superfluidity persists only if 



K>Kl 



TinL 



(96) 
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